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1.0.  BACKGROUND  AND  INTRODUCTION 
1.1.  Background 

Airborne  lasers:,  througn  measurement  ot  the  return  time  from 
the  ocean  surface,  ; f cord  a  method  tor  the  study  of  ocean  waves 
and  disturbances.  Theoretical ly,  given  a  sufficient  number  of 
irtasurements  snrpad  ’p  ,hi>-  -■  sr  the  ocean  Surface  u» Ci.  a  lu« 
flying  aircraft,  it  should  be  pcssiole  to  obtain  a  description  of 
the  surface  at  least  as  accurate  as  the  laser  measurer  enf . 
However,  trie  limited  freauency  of  measurements,  together  with  the 
fact  that  the  surface  is  continuously  chancing,  makes  toe  talcing 
of  a  "snapshot of  the  ocean  surface  exceedingly  difficult. 

Some  years  ago  the  U.  S.  Navy  sponsored  the  deve 1  cement  of 
an  airborne  laser  survey.  This  Navy  sponsored  won  .  References  5 
to  B,  generated  an  algorithm  ana  a  FORTRAN  computer  program  that 
took  time  parameters,  navigational  data,  aircraft  attitude, 
aircraft  altitude,  laser  beam  direction,  bottom  and  surface 
returns  to  generate  a.  representation  the  ocean  floor.  The 
underlying  concept  was  Optimal  < Ka  1  mar.  liter  Theory  which 

rsau.1  red  phrasing  this  estimation  problem  m  a  canonical  form  or 
model.  As  cart  of  this  model  was  a  simplified  representation  of 
wave  heights  as  Markoff lan  (exponentially  correlated)  noise. 

Thus,  implicit  in  tne  algorithm  was  an  estimate  for  wave  height 
However,  this  estimate,  reflecting  economies  in  model  ling,  could 
not  reflect-  the  full  potential  accuracv  inherent  in  the  laser 
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1.2.  Introduction 

A  contract  was  awarded  to  Papin  systems  AssonatPE .  Inc.. 

9  Seotetoer  l  c>36 .  for  tne  exploration  c?  tne  cctertiel  of  an 
airborne  1  aser  for  the  estimati  on  of  wave  ne  ic.hr.  The  aoprotcn 
teten  was  to  reproduce  old  results*,  :.s  onef  ly  described  in 
Section  2.  'It  soon  became  apparent  treat  the  simple  ana 
econom:  cal  three  or  four  state  models  used  .in  sari  ;er  wori  were 
i n a a equate  and  a  full  fifteen  state  model  for  a  P3~  aircraft  was 
emp  1 oyer  as  aescribed  m  Section  3.  Section  A  describes  the 
results  Detained  using  this  mocei .  The  results  of  vary mo  wind 
oust  stanc&rd  deviation,  wind  oust  correlation  time,  and  the 
stanca.ro  deviation  of  combined  wave  heipht  ar.a  measurement  error 
are  described  ,  in  Section  5.  '  As:-  .long  wave  lengths,  character  i  st  i 
or  deeo  ocean,  were  also  of  interest,  a  model  variation  was 
emo loved  to  separate  wave  height  anc  measurement  error,  as 
described  in  Sec  t :  on  t>.  The  effect  of  varying  pulse  repet  3  1 1  on 
rate  and  beam  rotation  rate  are  discussed  in  Section  3. 

Cone  1  us  ions  and  suggestions  for  possible  further  work  are  given 
i n  Sec  1 1 on  S.  -  ■  •  -  "  / 

2.0.  REPRODUCING  OLD  RESULTS 

War  I-  ceosn  by  reproducing  the  result  of  Reference  5.  The 
four  state  Attalt  Filter  Program  was  set  up  in  BASIC  in  an  IBM  X 
environment.  A  simoiified  measurement  simulation  was  set  uo 
i  ncorcor-a  t :  no  simulated  measurements  for  a.  2  degree  bias  error  i 
both  pitch  and  rc  11.  Aside  from:  this  pitch  and  roll,  there  were 
no  additional  simulated  errors.  (This  simple  simulation  is  nsec 
throughout  the  seauel .  .*  This  program  is  ref  erred  to  as  ATTALT2. 
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The  program  was  e  le-'C  1  sea  to  verify  correct  operation 
consistent  with  the  results  of  Reference  6.  It  soon  became 
evident  that  the  tour-  state  model  was  inadequate  for  two 
r  ease ns : - 

'.a’  If  a  bias  error  existed  in  pitch,  the  ATTALT2  Filter 
Program  would  comte  close  to  estimatino  this  Dias  when 
the  laser  is  fore  or  aft.  However,  as  soon  as  the 
beam  came  atnwar  t-sh  i  p .  the  estimate  would  drift 
toward  zero  consistent  with  the  hart  of f i an  model  for 
Pitch  error.  This  would  cause  unacceptable  errors  for 
wave  estimation.  A  similar  argument  pertains  to  roll. 
Thus,  at  least  two  additional  states  are  necessary  for 
D 1 1 c  h  a n d  roll  bias  e  r  r  o  r . 

(b )  I n  exerc i s i no  ATTALT2,  it  was  necessary  to  reduce  the 
standard  deviation  of  the  Ditch  and  roll  motion  to  0.2 
of  a  degree  and  the  vertical  motion  model  standard 
deviation  of  vertical  motion  to  0.001  meter  in  order 
to  achieve  accuracy  of  estimating  wave  height 
measurement  noise  to  the  c-der  of  0.1  meter.  Clearly, 
this  model  of  aircraft  motion  is  inconsistent  with 
real i tv. 

These  observations  led  to  the  conclusion  that  a  more  complex 
model  for  aircraft  motions  was  necessary.  After  some  flirting 
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3.0.  THE  LOCKHEED  P3C  AIRCRAFT  MODEL 

The  -first  foray  into  the  realm  of  full  aircraft  motion 
models  employed  equations  of  motion  from  Reference  2.  pp . 171-174 . 
There  was  concern  that  this  model,  for  a  jet  aircraft,  would  not 
oe  representative  of  the  P3C  aircraft  used  for  field  tests.  (As 
it  turned  out.  partial  results  with  this  model  did  not  differ 
significantly  with  results  using  the  P3C  model.)  Consequently , 
Reference  3,  a  set  of  differential  equations  for  the  P3C  aircraft 
was  obtained  through  the  good  offices  of  Mr.  John  W.  Clark,  Jr., 
Aerospace  Engineer  at  the  Naval  Air  Development  Center. 
Warminster,  PA. 

These  differential  equations,  together  with  their  phrasing 
in  the  canonical  form  of  Optimal  (Kalman)  Filter  Theory,  are 
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The  parameters  given  in  Appendix  B  are  referred  to  as  the 
baseline  model.  This  baseline  model  assumes  a  standard  deviation 
of  wave  height  plus  measurement  noise  of  0.3  meter  with  each 
measurement  having  independent  errors.  The  wind  gust  model 
assumes  a  standard  deviation  of  wind  gust  of  0.02  radian  (0.12 
degree),  a  moderate  condition.  Modifications  to  this  baseline 
model  were  made  as  described  in  the  sequel. 
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4.0.  RESULTS  USING  BASELINE  P3C  MODEL  ( ATTALT8) 

The  p-3C  Lacelme  model  was  analyzed  using  the  ATTALTB 
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A  plotting  program  was  developed  to  facilitate  plotting 
results.  The  following  discusses  these  plots  for  the  easel  me 
model „ 

4.1.  Plot  Of  Sigma  Of  Wave  Height  +  Noise  Estimate  For  2  Sec. 

Figure  1  is  a  plot  of  the  standard  deviation  of  the  estimate 
of  wave  height  plus  measurement  noise  as  they  are  lumped  together 
m  the  easel  me  model.  The  a  priori  standard  deviation  plus 
measurement  noise  is  0. 3  meter.  The  run  extended  for  two  sec on os 
and  points  were  plotted  after  every  fifth  laser  range 
me  a.  sure  mien  t . 

Note  that  the  standard  deviation  of  the  estimate  of  wave 
height  plus  measurement  noise  is  better  than  0.1  meter- - better- 
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without  knowing  wmch  of  the  two  this  sum  represents.  To  obtain 
an  e  ooro:.  l.mate  standard  deviation  of  the  wave  estimate,  the 
standard  deviation  of  wave  plus  measurement  noise  must  be  RMS " d 
with  the  measurement  noise.  That  is.  if  the  standard  deviation 
of  the  estimate  of  wave  height  plus  measurement  noise  is  0.08 
meter  and  that  of  the  measurement  error  is  0.1  meter,  then  the 
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resulting  standard  Ce.tatior.  teat  .  ave  height  is  r,.p  r  or  i  mat  e  J  > 

0.  128  me  ter  .  F1  .case  note  ttot  t  ne  atC'ce  er  simple  ?.  ino  i  ;es  toe  v, 
the  limiting  accuracy  ,  so  measurement  Error  went  to.  rero.  m 
estimating  the  wave  height  of  s  0 .  Z  meter  wave  would  he  0.12 
meter . 

The  plot  in  Figure  1  d  i  splays  considerate  le  fundamental 
frequency  rj th  s  0.2  second  period.  Part  of  tn is  can  he 
explained  by  the  two  degree  bias  inserted  in  pitch  arc  roll. 
However ,  removing  this  Pies  did  not  eliminate  the  fundamental 
component.  It  turns  out  that  the  mechanism  for  introducing  this 
fundamental  component  involves  a.  coupling  between  pi  ten  error  and. 
altitude  error-.  The  pitch  and  »]  t  etude  effects  are  acd  i  1 1  ve  when 
tne  laser  points  forward  but  subtract  wnen  pointed  aft.  Tne 
cancel  at  i  o  n  c  a  u  see  1  s.  r  g  e  ►*  u  r,  cert  a. :  n  t  y  w  hen  t  n  e  laser  p.  -  i  n  t  s  aft 
and  mini  mum  when  pointed  forward. 

4.2.  Plots  Of  The  Estimates  Of  Pitch  And  Roll 

F i cures  2  and  3  are  plots  of  tne  estimates  of  pitch  and 

trim  condition  t  and  e  and  an  es  1 1  rr.ate  of  "Mss"  which  wo*.  lc'. 
include  ecu l pment  mi saJ i on  merit  e^ror  s 1  us  in  the  case  of  pi  ten 
the  .  a  i  a  1 1  on  of  trim  condition  wi  tn  speed  and  load.  Note  that 
both  pitch  and  -‘oil  bias  estimates  aporoacn  2  degrees  <0.035 
radian;.  It  tales  considerable  time,  particularly  in  roll,  to 
cGiTie  up  with  the  estimate  of  bias.  However,  in  eacn  case,  the 
sum  cf  the  estimates  reoidl*  apo-oaches  two  degrees  reflect!  g 
the  fact  that  tne  accuracy  of  the  sum.  of  the  estimates  is  much 
better  than  that  of  the  components. 
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4.3.  Plot  Of  Sigma  Of  Wave  Height  +  Noise  Estimates  For  10  Sec. 

F  :  care  4  ceser.  ts  a-  clot  for  the  s  a  me  be  eel  me  co'.c?  i  t  a  r.r.=  a 
for  Fiiure  i  but  e  tenaed  to  ten  eeconpn.  F  ronts  a-*e  tic*’  ed 


t-  -8r'.  .  ■  '.'Z  as.,  cc-c  or  6-tj"  ■  ’ cecj  tre=  u  1  1'i't  at  i  or*  c  f  the  bean., 
i  t  seem 2  1 1 t e  two  different  curves  a'  e  Diet  ted  here.  It  on  1  , 


w  e  aits  t  h  ;7  t  W  a  .  •  1  ‘"I  *_■  .*  Cwt!"  C  ;J  !’  V  e  1  I  C*  Of*  j  :  t  Cl  f  h  e  I  r*  SEI 

C'Oint  i".c  Icwtr’i  which,  as  described  auo  re.  yields  a  i  nwer 

Tne  last  oemt  or.  this  curve  has  the  ta.Iut  o“  C .  OS  I  mete"  tc  nr. 

s  t  a  n  d  a  d  ae-'istion  of  the  est  imatc  c  f  wave  “.eight  :lu=  noise  .  i 

the  results  oescr  iced  below,  where  a  nous  parameters,  are  ,;>'iec 
it  i  =  this  last  point  at  the  end  ct  ten  seconds  that  is  c- lotted. 
5.0.  THE  EFFECTS  OF  VARYING  PARAMETERS 

The  parameters  of  the  wind  gust  me  cel  arc  the  magr  itu.de  of 
the  we  /65  clue  measurement  noise  were  -  aried  to  determine  the:*- 
effect  ;n  the  accuracy  of  wave  height  estimation. 

5.1.  Varying  The  Standard  Deviation  Of  Wind  Gusts  (ATTALT7X) 

A  series  of  runs  were  made  for  different  values  c-f  the 
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5.2.  Effect  Of  Varying  Wind  Gust  Correlation  Time  (ATTALT7Y) 
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5.3.  Varying  Wave  Height  Plus  Measurement  Noise  (ATTALT7Z) 
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6.0.  SIMPLIFIED  MODEL  FOR  LONGER  WAVES  ( ATTALT9) 

the  p  a  e  e  1  1  n e  ~  jce-1  :  OEi:uiT|EE  ot  es cr  ]  a* a e '  ranee 
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■ho  baseline  mooe  I ,  the  dirts,  roe  on  the  euc  sce  between 


<T>ee ~ u remen t s  is  appro  lmetely  2c  meters. 


or  sumo t  :  or  i  r.  cep  eocene  e  is  a  valid  one  for  shallow  seas . 
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iT'O  t  e  s .  "  r  e  o' ore  the  effect  cf  wave  length  on  wave  height 
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7.0.  VARYING  PULSE  REPETITION  RATE  AND  VARYING  BEAM  ROTATION 
RATE 

three  r.jr.c  were  rosde  tc  explore  the  lmc-act  on  a::ur-.c;.  oi 


mo  1  as 


'he  first  :  i  r 


r  ep  etit  ion  <-ate  and  beam  rote  t  ion  rate, 
a  <r DC  ■  *  1  ca  1 1  ort  of  ATTALT8  to  increase  the 


Dulse  reset: ton  rate  from  4 00  to  600  suises  oer  second .  This 
cnarpe  reduced  the  standard  deviation  c:  the  error  in  est  ’.mat  ion 
of  *4*  v-e  r.s  iciht  -*•  measurement  error  from  1.0022  tc  0 . 0759  meter. 

The  second  run  was  a  modi  float  ion  of  ATTALT8  to  increase 
beam  •■ctetion  rate  from  five  to  ten  revolutions  per  second.  This, 
chan  do  educec  the  standard  deviation  of  the  estimate  of  wave 
he:  g  fit  plus  tne  error  in  measurement  from  0.0G22  to  0 . 0802 
meter  .. 

The  third  run  incorporated  Doth  of  these  chances 


s  i  mu  t  an  ecus  .1 


i e  pl;. i  se  r  tr1  o a  f  i  1 1  on  rate  was  increased  t c  vu 'J  j 


pulses  per  second  and  the  beam  rotation  rate  was  increased  tc  ter. 
revolutions  per  second.  This  cnanoe  reduced  tne  standard 
deviation  of  the  estimate  of  wave  height  plus  measurement  noise 
f r o m  0 . 0 G 22  t o  •  1 . C "3o  m e t e r . 

It  should  be  remembered  that  in  order  to  obtain  the  standard 


oev l s t i on 


RMS  '  c  wi 


h e :  oh 


the  estimate  of  wave  height,  these  number'  must  be 
re  measurement  error  standard  deviation  cf  0.1 
the  standard  deviation  of  the  estimate  c~  wave 
e  above  three  runs  is,  respectively:  0. 125,  0. 1 28, 


and  0. 124  meter.  This  is  not  much  different  than  tne  baseline 
results,  cf  0.  12-  meter.  Thus.  as  long  as  measurement  error 
remains  at  0.  1  meter,  it  is  tne  control  ling  source  of  error  and 
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mcreasec  Dulae  rate  and  beam  rotation  rate  does  not 
significantly  improve  the  estimate  of  wave  height. 

8.0.  CONCLUSIONS  AND  RECOMMENDATIONS 

For  baseline  conditions,  where  the  laser  is  the  prime  roll 
and  pitch  reference,  a  shallow  water  environment,  and  the  wave 
height  has  a  standard  deviation  of  about  0.3  meter,  it  is  clear 
that  wave  height  can  be  determined  to  the  order  of  the  accuracy 
of  the  measurement,  0.1  meter.  Reducing  the  measurement  error 
will  somewhat  improve  the  accuracy  of  the  estimate  of  wave  height 
under  these  c ond i t i ons . 

In  the  shallow  water  environment,  accuracy  deteriorates 
rapidly  for  larger  waves.  Reducing  the  measurement  error  would 
have  little  effect-  under  these  conditions. 

In  the  deep  ocean  environment,  with  longer  wave  lengths,  the 
accuracy  of  wave  height  estimation  decreases  significantly.  It 
should  be  pointed  out  that,  m  this  long  wave  case,  the  errors  in 
successive  estimates  would  be  correlated  yielding  a  good  picture 
of  tne  shape  of  a.  wave  but  not  its  absolute  magnitude. 

If  the  above  accuracies  predicted  in  this  report  are  not 
adequate,  if  might  be  necessary  to  explore  the  introduction  of  a 
stable  reference.  This  might  take  the  form  of  an  autopilot,  an 
Attitude  and  Heading  Reference  System,  or  an  accurate  inertial 
system.  For  the  baseline  condition,  significant  improvement 
would  only  be  obtained  by  also  reducing  measurement  error. 

It  is  recommended  that  the  several  avenues  for  improved  wave 
height  estimation  be  further  studied  in  the  light  of  Navy 
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FIGURES  USED  IN  SECTIONS  1  TO  7 


FAGIN  SYSTEMS  ASSOCIATES,  INC.  -  FIGURE  1.  PLOT  OF  S.D.  OF  lil+MEAS  ERROR  ESTIMATE 
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FAGIN  SYSTEMS  ASSOCIATES,  INC.  -  FIGURE  3.  PLOTS  OF  ESTIMATES  OF  PHI  &  PHIB 
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Laser  Wave  Estimation 

NSTL  Contract  N00014-88-C-603: 

Prepared  by:  F3A  Inc.,  Wynnewooc.  PA 


inininwmininTiinininsnininTiinininwu'Hn 
r-  -  oj  r-  cu.  r'-  .  cu  cu  -r^-ioj  -r^ico 

.  CU  -  -rH  -  -  -  I  -  •  OJ 

CU  CU  -H  »  1  -HIt-1  cu  i 

i  i  i 

T-30T  x  < - S3313U  -  o»oyuois 


o 

UJ 

O') 

e 

LlI 

<5> 


yc 


a 

i 


LU 

2 Z 


-IS- 


-2.75 


Laser  Wave  Estimation 

NSTL  Contract  N00014-BS-C-6031 

Prepared  Dy :  rSA  Inc.,  Wynnewood,  Pri 


REFERENCES 


Fa'.,  i  n,  S.  I  .  (  j  ?w‘'i  1  .  "Recure  i  ye  Linear  F.'eareesi  on 

7  he  or  *  .  Dpt  3  mal  Filler  Theory  ,  an  Error  Analyses  of 
Opt  m.a  I  Evstemr.  "  IEEE  I  n  ternat  i  ona  1  Conv,  Rec  .  , 

Vo  1 .  '  L.  F't .  1,  po. 

F'r\  per. ,  E.  D.  ,  Giant  ord  IJniv.  ,  ana  Ho,  Y .  -C .  .  Harvard 
ljn  i  v  .  ,  '1  v/,o  >  .  '•  App  l  i  ed  Up  1 1  ma  3  Con  t  ro  1  , 

Dp  t  i  ni  i  z  a  1 1  on  .  Fet  i  ma  1 i  on  a  n  d  C  on  trol  .  "  G  i  n  n  a  n  d 
Comp an,  .  pp  17]  —174. 

Lear  Si  pier,  Ir.c  .  .  As  t  ron  i  c  s  Div.  i  17*6°)  .  "  AFCS 
S  i  mi1  i  a  t  i  on  Report  For  Lockheed  P-3C  Aircraft.  " 

Prepared  for:  Loc  Dheed- Ca  1 1  torn  l  a. .  Burbank,  CA. 

Da  la  Item  No.  05-241-00. 

A vc a  Everett  Research  Laboratory ,  Inc.  < 19753 . 

"Aji'Scrne  Oceenogra.pl' •  lc  L I DAR  System."  Prepared  for: 
NASA,  Wallops  flight  Center.  Contract  No.  NA56-2653. 

Fag i n ,  S.  L. .  F  a  a i n  Assoc i eves.  (June  1 c77) .  "Model  for 
Op 1 1 me  1  Laser  Curve  -mg."  Prepared  for:  Naval 

uceonog  rapn  i  c  Of  1  ice.  Contract  No.  N684p3-7c-C-007. 

Byrnes,  H.  J.  ,  MC*<G  Development  Group,  Naval 
Ooeanogreph  i  o  L.  ah  or  a  t  ory,  and  r:  a  gin,  5.  L.  ,  Fa  gin 
Associates.  (lc73>.  "Optimal  Filtering  and  Analysis 
—  f  oc  enn  i  no  uaser  Data.  "  NOkDA  Technical  Note  24, 

June  1978. 

Fag  in.  S.  L. ,  Fag  in  Associates.  U.979) .  "HALS 
Shipboard  Curve,  Processor .  "  Preoared  for':  Naval 
Dc eanog r aph i c  Office.  Contract  No.  N6B463-7S-C-0041 , 
24  Ap  i’il  1979. 

Byrnes.  H.  J.,  MCLG  Development  Group,  Naval 
Oceanographic  Laboratory.  (1982).  "HALS  Post 
Processing  Software  Design."  NORDA  Technical  Note 
161,  NSTL  Station,  MS,  September,  1982. 


Laser  have  Estimation 

NSTl  Contract  N-.»001 4-85-C-6031 

Frepareo  by:  F5A  Inc.,  h/nnewood,  F’A 


APPENDIX  A. 

THE  CANONICAL  MODEL 

in  order  to  employ  the  theory  and  aloon  trims  of  Optimal 
u. a  i  man)  Filter-  Theory,  it-  is  necessary  to  place  a  problem  in  a 
tDf?c  i  a  1  form  called  "The  Canonical  Model".  A  brief  description 
of  this  mode i  is  given  here  for'  completeness. 

In  the  sequel,  the  notation  of  Reference  1  is  adopted. 
L"iefly,  the  error  model  is  assumed  to  De  represented  by  the? 
first  order  differential  eauations: 

X  =  A(t)X  +  N(t)  ( A  —  1 > 

where: 

X  is  a  "k"  dimensional  column  "state"  vector. 

A(t;  is  a  I  ::  k  "system  ’  matri;:,  which  may  vary 
with  time. 

M(t)  is  a  " k "  dimensional  column  vector  of  white  noise 
inputs  such  that 

E  -C  IM  ( t  ^  )  N  f  (  1 2  )  i  =  n  £  ( 1 1  - 1-  2  >  Q  <,  1 1  1 
and  E  C N  <  t  >  j  =  0 

where  £  is  the  Dirac  deltai  function 

q(tj;  is  a  k  k  "noise"  matri;:. 

In  the  modelling  adopted  here,  the  state  vector  X  includes 
"noise"  states  sucn  as  measurement  errors  and  error  sources. 

Most  frequently  q  is  constant  and  diagonal,  the  diagonal  elements 
being  tne  single-ended  power  spectral  density  of  the  white  noise 
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necessary  to  create  the  no : :e  state  of  the  dee: red  variance. 

At  times  tx ,  "r"  dimensional  observations  Ytt, >  are  tel  en 

defined  h v  tne  eauaticn 

Y  <tx  ;  =  M  ( t  j  )  X  t  i  )  ( A- 2  ) 

M  .  t )  is  a  measurement  of  dimens  . on  r  ,  l  where  r  I-  .  M(t,  ; 

i.  j. 

not  necessarily  constant  but  is  assumed  to  piecewise  eontmuou 
Equations  A— 1  and  A-2  can  be  represented  by  the 
multidimensional.  Canonical  E-iocl  diaoram  of  figure  A-l .  In  t 
figure,  double  lines  indicate  the  flow  of  vector  quantities,  a 
niotfi ;;  bores  multiply  their  vector  inputs  from  the  left. 

The  optimal  filter  or  smoothing  problem  is  to  obtain  a 
"best"  estimate  for  the  state  vector  " K "  bases  on  these 
o&ser  rations  and  any  prior  knowledge  of  the  states  summed  up  h 
a  n  i  n  i  1 1  a  1  c  o  v  a  r  i  a  nee  m  a  tr  i P  ( 0  -  >  . 

To  specif-  a.  model,  it  is-  vherefor  necessary  to  specify  A 
a,  P  ( 1  .  •  ( 0-  •  ,  an:  tne  measurement  ma.tr  ices. 


N(f)s=C> 
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Figure  A-l 
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APPENDIX  B. 
THE  P3C  MODEL 


'i  he  differential  equations  for  the  P3C  aircraft  motion  about 
trie  cond:  tion  m  straight  and  level  flight  wart  obtained  from 
Reference  3..  Flight  Condition  1C  was  chosen  as  most  closely 
represent  mg  the  situation  of  interest.  Standard  aerodynamic 
nomenclature  is  employee).  The  reader  is  referred  to  the 
reference  for  definition  of  terms.  No  effects  of  control 
surfaces  were  included:  however,  additional  states  were  added  as 
shown  in  Figure  B-l,  as  follows: 

pg  =  bias  in  Ditch  either  due  to  the  trim  of 
the  aircraft  and/or  misalignment. 

op  =  bias  in  roll  either  due  to  trim  of  the 
aircraft  and/or  misalignment. 

w  =  w a v e  h e i o h t  p 1 u s  m e a surement  e r r o r 
modelled  as  a  Marrcftian. 

In  addition,  the  variables  aoust  and  JiQUS-  were  modelled  as 
Marrof f i an  random  variables.  Tne  resulting  linear  differential 
equation  of  dimension  15  is  shown  in  Figure  B-l. 


he 


shown  in 


correspond i ng  state  vector  "X"  and  the 
Figure  B-2.  The  "q"  matrix  is  defined 


"A" 
by : 


matrix  are 


q  (5,5)  —  u.  i  )  /  71  TQUEt 


q (12, 12 i  =  2  m^<oQUst) 


/  TT  T 
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q  (  1  4  ,  1  4  )  =  2  <7  *-  (  w  >  /  Tt  t  w 

Air  other  values  of  the  matri::  are  zero. 

The  initial  covariance  ,natr  l is  diagonal  with  these  having 
standard  deviations  respectively  of  10,  0.1.  0.1,  0 .01 ,  0. 02, 

0 . 1 .  0.1,  0 . 0 1 ,  0 . 0 1 ,  0 . 0 1 ,  0 -  1 ,  0 . 02 ,  0 .  1 ,  0 . 3 ,  0 . 5 . 

A  list  of  the  reme i n i no  parameters  for  the  baseline  P3C 
model  is  given  in  Figure  E--3.  All  units  are  Meters,  Radians  arc 
Sec  on d s . 
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Figure  B-l.  P3C  DIFFERENTIAL  EQUATIONS 
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APPENDIX  C. 

RELATIONSHIPS  BETWEEN  ERRORS  IN  VERTICAL, 
ALTITUDE  AND  WAVE  HEIGHT 


A  :  r. i  :  t  itu  c  :  so raiTi ■  s h u w 3. n c-  t  n  e  rs  1  a 1 1  on s* h  i  p> ~  d e t w tr e n 
aititj.CE-  »H)  ,  the  distance  the  laser  beam  t  rave  is  to  the  surface 
;P  -  .  and  the  an  ole  (A>  between  the  beam  and  the  vertical ,  is. 
h a w n  in  P" i 0 u r e 
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It  follows  that  if  wave  heioht  error  is  0.1  Meter,  then  dH 
plus  134  dA  must  add  up  to  0.1  lister .  Unless  tney  are  highly 
correlated,  this  implies  that  dH  is  the  order  of  0.1  Meter  and 
dA  is  the  order  cf  2 .  t  minutes  of  arc. 
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